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摘 要 : 为 提高 生长 初期 低 覆 盖 度 作物 长 势 的 还 感 监测 精度 ， 需 要 消除 灌溉 引起 的 土壤 水 分 背景 变化 对 归 一 
化 差 值 植被 指数 (NDVD 的 影响 。 为 了 实现 棉花 生长 初期 灌溉 信息 提取 与 校正 ,提高 棉花 作物 长 势 监测 与 产量 
预 判 精度 ， 本 文 以 美国 加 利 福 尼 亚 州 San Joaquin Valley 的 2 个 棉花 地 块 为 研究 区 ， 选取 棉花 生长 初期 灌溉 过 
程 中 的 遥感 影像 ， 构 建 两 种 灌溉 信息 提取 方法 (分 阶段 国 值 法 和 灌溉 线 提取 法 )， 确 定 最 优 灌溉 像 元 提取 方法 ; 
比较 分 析 灌溉 与 未 灌溉 情况 下 棉花 的 NDVI 与 归 一 化 差 值 水 分 指数 (NDWD 以 及 土壤 调节 植被 指数 的 关系 ， 提 
取 含 有 灌溉 信息 的 像 元 ,并 对 NDVI 进 行 校正 ， 消 除 灌 溉 对 NDVI 的 影响 .研究 结果 表明 : 在 棉花 生长 初期 ， 灌 
溉 与 未 灌溉 像 元 NDVI 变化 率 达 12%， 差 异 较 显著 ; 灌溉 与 否 的 棉花 NDVI 与 NDWI 间 均 存在 极 显著 的 线性 
关系 ,决定 系数 在 0.80 以 上 ; 利用 灌溉 线 方 法 提取 灌溉 信息 与 分 阶段 国 值 相 比 精度 更 高 ， 精 度 达 88% 以 上 ; 校 
正 后 线性 回归 模型 精度 达 0.95， 灌溉 校正 效果 明显 ， 灌 溉 与 未 灌溉 像 元 的 NDVI 差异 减 小 至 2%。 本 研究 通过 
对 含有 灌溉 信息 像 元 NDVI 值 的 校正 ， 去 除 灌溉 对 NDVI 造成 的 影响 ,反映 了 真实 的 植被 信息 ， 可 实现 对 作物 
生长 初期 长 势 的 准确 遥感 监测 ,为 还 感 定量 监测 提供 便利 。 

关键 词 : 棉花 ; 生长 初期 ; 灌溉 信息 ; 植被 指数 ; 归 一 化 差 值 植被 指数 (NDVI); 归 一 化 差 值 水 分 指数 (NDWD 
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Using remote sensing to extract and correct irrigation data 
during early cotton growth stage 
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Abstract: Vegetation index ls affected by background soil, especially in the early stage of crop growth. When vegetation cover is low, 
the effect of background soil is very obvious. In order to improve the precision of remote sensing (RS) monitoring on crop growth in the 
early growth stage, it ls necessary to eliminate the effect of background soil moisture due to irrigation on normalized difference vegeta- 
tion index (NDVD. Agricultural irrigation districts have failed to develop an effective method to eliminate difference in NDVI change, 
which has in turn hindered efforts to limit the effect of irrigation on NDVI. Thus, in order to increase the accuracy of RS monitoring 
of crop growth at early stage, this study explored the effects of difference in soil moisture information between irrigated and 
non-irrigated cotton field on NDVI. Two cotton plots in San Joaquin Valley in California (US) were selected as the research area. Day 
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of Year (DOY) 174 was determined as the critical phase at early growth stage of cotton for the extract of irrigation data through band 
reflectance, NDVI analysis of cotton field for 2002. Based on RS images, NDVI, normalized difference water index (NDWD, soil 
adjusted vegetation index (SAVI) and modified soil adjusted vegetation index (MSAVD of irrigated and non-irrigated pixels were 


calculated. Also the relationships between NDWI and different vegetation indexes (VIs) were analyzed, and the two methods [the 


standard deviation of the NDWI method (STDwn and irrigation line extraction method (based on relationship between NDVI and 


NDWI of irrigation and non-irrigation pixels, IR_L)] were used to extract the irrigation data. Then the accuracies of different meth- 


Ods were compared to determine the optimum extraction method of irrigation information. The IR_L method was next used to extract 


irrigation data and correct the NDVI of irrigation pixels in the early stage of cotton to improve monitoring accuracy of cotton growth. 


The results showed that difference in NDVI between irrigation and non-irrigation pixels was as high as 12% in the early growth of 


cotton. There was an extremely significant linear correlation between NDVI and NDWI of both irrigation and non-irrigation pixels, 


with coefficients of determination greater than 0.80. Compared with STDwr method, IR_L method had a higher accuracy and with a 


precision greater than 88%. Through IR_L model correction, the accuracy of irrigation linear regression model was as high as 0.95. 


With this, correction effect of irrigation was obvious and the difference in NDVI between irrigated and non-irrigated pixels dropped 


to 2%. Thus in this study, NDVI with irrigation data was corrected, the effect of irrigation on NDVI eliminated while the effect of 


background soil moisture reduced. Finally, the study reflected the true vegetation data, obtained accurate remote sensing monitoring 


of cotton growth at the early growth stage and provided convenient monitoring method of crop growth via remote sensing. Moreover, 


it promoted accurate irrigation towards saving water resources. 


Keywords: Cotton; Early growth stage; Irrigation data; Vegetation index; Normalized difference vegetation index (NDVD); 


Normalized difference water index (NDWI) 


遥感 技术 具有 宏观 、 客 观 、 快 速 、 低 成 本 等 特点 ， 
成 为 近 几 十 年 作物 估 产 和 长 势 监测 的 新 兴 技术 5 J。 及 
时 准确 地 监测 作物 生长 初期 植被 长 势 ， 可 以 对 作物 
产量 信息 进行 预 判 8B4。 目 前 ,在 农业 遥感 监测 中 ， 
植被 指数 (VD 应 用 广泛 694， 国内 外 用 于 长 势 监测 的 
植被 指数 有 多 种 ， 其 中 归 一 化 差分 植被 指数 (NDVD 
使 用 最 为 广泛 。 裴 志和 远 等 中 利用 同一 生长 季 不 同 生 
育 期 及 不 同 生 长 季 同 一 生育 期 的 NDVI 值 ， 构 建 不 
同 的 NDVI 时 间 序 列 曲线 的 特征 参数 ,建立 了 作物 
长 势 监测 模型 ; 黄 青 等 中 利用 MODIS_NDVI 分 析 确 
定 了 冬小麦 (Triticum aestivum L.) 信 息 提 取 的 NDVI 
阅 值 ， 建 立 了 冬小麦 面积 提取 模型 ， 得 到 2011 年 冬 
小 麦 整 个 生育 期 长 势 情况 。 

然而 ， 由 于 研究 区 的 气候 、 水文 条 件 不 同 , 土壤 
背景 会 存在 差异 。 植 被 指数 会 受到 农田 土壤 背景 影 
响 M' "1 尤其 是 在 作物 生长 初期 ,植被 覆盖 度 较 低 ， 
土壤 背景 的 影响 更 加 明显 。 土 壤 背 景 对 植被 造成 的 
影响 与 作物 对 地 覆盖 度 一 一 叶 面 积 指数 (LAT) 关 系 
密切 ,覆盖 度 越 大 ,土壤 背景 影响 越 小 。 因此 在 研究 
土壤 背景 对 NDVI 的 影响 时 ,不 能 忽略 作物 对 覆盖 
度 的 影响 。 刘 明 等 (1 根据 NDVI 对 叶 面 积 指数 进行 反 
演 ,， 验 证 NDVI 与 LAI 之 间 存 在 显著 相关 性 。 作 物 生 
长 初期 覆盖 度 较 低 , NDVI 值 小 于 0.5, 土壤 背景 影响 
较 明 显 ， 需 要 消除 土壤 背景 对 NDVI 的 影响 (1。 

Rondeaux 等 (站 从 NDVI 出 发 ， 对 比 了 土壤 调节 植 
被 指数 等 多 种 植被 指数 消除 土壤 背景 影响 的 能 力 ， 
指出 在 不 同 环境 条 件 下 ， 不 同 植被 指数 消除 土壤 背 


景 影响 的 能 力 不 同 。 然 而 土壤 调节 系列 的 植被 指数 
应 用 虽 广 泛 ,物理 意义 明确 ; 但 由 于 土壤 调节 参数 
确定 存在 一 定 的 困难 ,不 能 完全 消除 土壤 等 背景 干 
扰 。 对 于 灌溉 农业 区 ,在 整个 作物 生育 期 进行 灌溉 
处 理 ， 导 致 作物 生长 期 内 灌溉 与 未 灌溉 区 域 么 感 影 
像 呈现 不 同 特征 ,造成 NDVI 值 的 显著 差异 09。 因 为 
灌溉 造成 土壤 水 分 信息 分 布 不 均匀 ,土壤 水 分 差异 
会 造成 土壤 背景 亮 暗 的 程度 不 同 ,对 于 相同 的 冠 层 
结构 ， 暗 土壤 背景 的 冠 层 NDVI 值 大 于 亮 土 壤 背 景 
的 冠 层 NDVI 值 。 目 前 ， 针 对 水 分 遥感 监测 更 为 直接 
的 指标 是 归 一 化 差 值 水 分 指数 NDWD。 很 多 学 者 利 
用 NDWI 进 行 水 分 以 及 土壤 水 分 信息 的 提取 。 程 晓 娟 
等 0 指出， 利用 改进 NDWI 监 测 作物 水 分 含量 ， 对 
于 评估 作物 水 分 盈亏 平衡 、 指 导 农 业 生产 灌溉 具有 
重要 意义 ; Lu 等 9 也 通过 HJ-1A/B 影 像 提 取 NDWL， 
a a J 分 析 , 精度 较 高 。 然 而 目前 
研究 只 是 针对 土壤 水 分 进行 提取 以 及 对 作物 水 分 合 
量 进行 监控 ， 对 灌溉 农业 区 的 作物 长 势 监测 少 有 提 
出 有 效 消除 NDVI 差 异 变 化 的 方法 , 没有 从 根本 上 
解决 灌溉 对 于 NDVI 造 成 影响 的 问题 。 

为 提高 生长 初期 作物 长 势 遥感 监测 的 精度 ， 深 入 
研究 灌溉 与 未 灌溉 情况 下 土壤 背景 水 分 信息 差异 对 
NDVI 的 影响 ， 本 文 以 棉花 (Gossypium sp.) 为 研究 对 
象 进行 如 下 研究 : 1) 为 消除 灌溉 对 土壤 背景 的 影响 ,分 
析 棉 花生 长 初期 灌溉 与 正常 状态 下 的 各 种 植被 指数 的 
差异 ; 2) 利 用 Landsat 30 m 空间 分 辩 率 遥 汪 感 影像 ， 获取 
NDVI 及 反映 土壤 背景 的 植被 指数 时 间 序 列 ， 进行 定 
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量 分 析 ， 明 确 表征 灌溉 后 土壤 背景 水 分 状况 的 植被 指 
数 ; 3) 确 定 灌溉 信息 提取 方法 以 及 NDVI 的 校正 方法 。 
本 文 研 究 结果 对 于 低 履 盖 度 下 的 NDVI 应 用 具有 一 定 
参考 价值 ， 有 利于 促进 生长 初期 作物 长 势 精准 监测 ， 
为 遥感 定量 化 与 精准 农业 发 展 货 定 基础 。 
1 数据 与 方法 
1.1 ”研究 区 域 概况 

1 为 研究 区 位 置 及 影像 图 。 本 文 研究 区 的 两 
个 地 块 种 植 的 农作物 均 为 棉花 ,位 于 美国 加 利和 福 尼 
亚 州 南部 的 San Joaquin Valley 西 侧 (图 1a)， 属 河流 
三 角 洲 地 区 , 土壤 为 棕 壤 土 ， 整个 农场 面积 约 为 
7 000 hm2， 每 个 地 块 面积 约 为 60 hm*( 图 1b)。 该 区 
属于 地 中 海 气候 ， 夏 季 炎 热 和 干燥 ,冬季 湿 润 凉爽 。 雨 
季 主 要 集中 在 11 月 至 次 年 4 月 , 但 降雨 量 极 少 , 年 
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1 研究 区 位 置 及 影像 图 (a: 加 州 ; b: San Joaquin Valley 农场 ; c: A 地 块 174 d 假 彩 


平均 降水 量 约 70 mm。 光 热 资 源 丰富 ， 最 高 温度 达 
40 人 以 上 , 蒸发 量 大 ， 屋 夜 温 差 在 16 左右 ， 有 利 
于 棉花 的 生长 , 是 棉花 的 高 产 区 1"。 研 究 区 农场 在 
棉花 整个 生长 期 实施 灌溉 ， 在 灌溉 初期 为 喷灌 ， 此 
后 每 2~3 周 进行 沟 灌 。 图 lc 与 图 1d 为 A 地 块 2002 
年 第 174 d 与 190 d 假 彩色 影像 , 图 le 为 B 地 块 174d 
假 彩 色 影 像 。 图 1c/e 可 以 看 出 地 块 的 灌溉 方式 是 从 
左 至 右 、 从 上 至 下 灌溉 ， 灌 溉 与 否 差异 十 分 明显 ; 从 
1d 的 第 190 d 影像 看 出 ,地 块 无 明显 灌溉 差异 ， 
但 该 地 块 在 相同 区 域 , 174 d 时 存在 灌溉 差异 明显 界 
线 。 选 择 190 d 是 因为 此 时 研究 区 棉花 地 块 灌溉 比 
较 均 一 , 与 174 d 形 成 对 比 。174 d 与 190 d 相差 16 d， 
由 于 本 文选 择 的 遥感 影像 是 Landsat TM 5 和 
Landsat ETM 7 两 种 影像 构成 的 时 间 序 列 ， 时 间 分 
辨 率 为 16 d。 


影像 ; 


d: A 地 块 190 d 假 彩 色 影 像 ; e: B 地 块 174 d 假 彩色 影像 ) 
Fig.1 Maps of the location and images of the study area (a: California map; b: San Joaquin Valley Farm map; c: false color image of 
plot A on DOY 174; d: false color image of plot A on DOY 190; e: false color image of plot B on DOY 174) 


1.2 ”影像 获取 与 预 处 理 

本 研究 主要 研究 棉花 灌 波 时 期 遥感 影像 ， 根 据 
美国 棉花 生长 期 的 特点 ， 下 载 研究 区 2002 年 棉花 
梨 土 期 与 生长 期 的 Landsat 5 TM 和 
Landsat 7_ETM 的 时 间 序 列 遥 感 影像 ， 时 间 分 辨 率 
为 16 d,， 空间 分 辩 率 为 30 m。 对 影像 数据 通过 
ENVI5.1 进行 辐射 定 标 、 大 和 气 校正 ， 并 按照 研究 区 
范围 矢量 图 进行 裁剪 。 


1.3 ”植被 指数 选取 与 计算 

NDVI 是 植被 长 势 监测 应 用 最 广泛 的 植被 指数 。 
土壤 调节 植被 指数 (SAVD 是 NDVI 的 一 个 改进 , 能 
减 小 土壤 等 背景 信息 ,增加 植被 信号 。 调 整 型 土壤 
调节 植被 指数 (MSAVD 是 众多 土壤 调节 植被 指数 中 
最 简单 的 ， 避 免 了 土壤 线 和 土壤 调节 参数 的 获取 。 
NDWI 利用 近 红 外 与 中 红外 进行 波段 运算 ， 用 于 研 
究 植 被 的 含水 量 ""， 本 文选 取 NDWI 对 灌溉 信息 进 
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行 研 究 。 各 植被 指数 计算 公式 如 下 : 


Nv MR RED 上 
NIR+RED 
是 (+ 万 ) (2) 
NIR+RED+L 

2NIR +1— (QQNIR +1)* -8(NIR -RED 
MSAVI= V | ) (3) 
NDWI= MR-SWR (4) 

NIR + SWIR 


式 中 : RED、NIR、SWIR 为 红 波段 、 近 红外 、 短 波 
红外 反射 率 。 工 为 土壤 调节 参数 ,描述 了 土壤 反射 特 
性 ， 当 植被 覆盖 度 较 低 时 ,推荐 L 值 为 1.0; 当 植 被 
覆盖 率 中 等 时 ,推荐 L 值 为 0.5; 当 植 被 覆盖 度 较 高 


a 


0 
126 142 158 174 190 206 222 238 254 270 286 302 


积 日 DOY 


First derivative 


阶 导 数 


| -0.010 


1 
—0.005 


—0.015 


-0.020 


时 , 推荐 工 值 为 0.25。 
1.4 ”棉花 灌溉 信 息 研 究 时 相 的 选取 
大 田 作 物 冠 层 光 谱 反 射 率 以 及 土壤 背景 对 其 造 


成 的 影响 程度 与 作物 覆盖 度 ( 叶 面 积 指数 ) 大 小 相关 。 
覆盖 度 越 大 ,土壤 背景 影响 越 小 。 本 研究 主要 研究 


作物 生长 期 的 灌溉 影像 ， 通 过 时 间 序 列 曲线 的 特征 参 
数 得 知 (图 2a), 棉花 从 174 d NDVI 上 升 速 率 增 加 ， 
NDVI 值 低 于 0.5, 由 NDYVI 与 叶 面 积 指数 的 相关 性 得 
知 174 d 植被 覆盖 度 并 不 是 很 高 ， 特 选取 174 d(2002 
年 6 月 23 日 ) 为 研究 灌溉 信息 关键 时 期 , 既 考 虑 到 作 
物 对 地 覆盖 度 的 影响 ， 也 充分 显现 出 土壤 背景 对 
NDYVI 的 影响 。 


0.020 
0.015 
0.010 


0.005 


积 日 DOY 


图 2 棉花 生长 期 NDVI 时 间 序 列 曲线 (a) 及 NDVI 一 阶 导数 曲线 (b) 


Fig.2 NDYVI time series curve (a) and NDVI first derivative curve (b) during cotton growing period 


1.5 ”灌溉 像 元 提取 及 校正 方法 
1.5.1 标准 差 方 法 

标准 差 能 反映 一 个 数据 集 的 离散 程度 ， 未 灌溉 水 平 
下 的 标准 差 反 映 了 正常 状态 下 作物 NDWI 的 离散 程度 ， 
本 文选 取 未 灌溉 像 元 NDWI 的 标准 偏差 (STDEVNwwi) 作 


为 分 阶段 阅 值 进行 灌溉 信息 提取 ,依照 公式 
dn wi = es 分 别 根据 NDWI+ dnpwi 进行 分 


梯度 NDVI 灌溉 信息 提取 , 简称 为 STDw 方 法 。 
1.5.2 ”灌溉 线 法 

采用 灌溉 线 进 行 灌溉 像 元 提取 ， 简 称 为 IR 工 方 
法 。 由 NDVI 与 NDWI 之 间 的 散 点 图 (图 3) 可 以 发 现 ， 
灌溉 与 未 灌溉 的 散 点 图 之 间 存 在 明显 的 分 界线 ， 本 
文 将 其 定义 为 灌溉 线 。 在 覆盖 度 逐 渐 增 加 时 ， 即 
NDVI 增 大 ,未 灌溉 像 元 的 趋势 线 与 灌溉 像 元 的 趋 
势 线 之 间 的 距离 会 逐渐 变 小 ， 灌溉 线 方程 的 斜率 采 
用 灌溉 趋势 线 与 未 灌溉 趋势 线 的 平均 值 ， 截 距 采 用 
两 条 直线 的 截 距 的 平均 值 ， 可 以 在 一 定 程度 上 降低 
高 植被 覆盖 度 下 土壤 背景 信息 弱化 的 影响 。 


NDVI= nr x NDWI+ nr (5) 


式 中 : Ki、Kui 分别 为 灌溉 、 i NDVI 与 
NDWI 线性 模型 的 斜率 ，b;.、bni 为 灌溉 像 元 、 未 灌溉 
像 元 散 点 图 的 所 在 直线 截 距 。 
1.5.3 ”灌溉 校正 方法 

本 文 的 目的 是 消除 灌溉 对 NDVI 的 影响 ,将 
灌溉 像 元 的 NDVI 调 整 到 未 灌溉 状态 的 NDVI, 使 
灌溉 前 后 植被 指数 具有 可 上 比 性 。 因 此 ,根据 调整 
原理 提出 校正 方法 , 分 别 以 灌溉 与 未 灌溉 各 像 元 
趋势 线 的 斜率 与 截 距 为 参数 进行 校正 ， 简 称 为 
IR LL 方法 。 


NDVIR r = 


天 ji 十 天 jj bi. 十 Di 
2 


(6) 


式 中 : NDVIrL 为 调整 后 的 未 灌溉 棉花 NDVL 
NDVIm 为 灌溉 后 的 NDVI，Ki、Kwir 分 别 为 灌溉 、 

未 灌溉 状态 下 的 NDVI 与 NDWI 线性 模型 的 斜率 ， 
pr、o 为 灌溉 像 元 、 未 灌溉 像 元 线性 模型 的 截 距 。 


天 
人 入 NDVYIR 十 (Or biir) 
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D 灌溉 像 元 Irrigated pixels 
A ”未 灌溉 像 元 Non-irrigated pixes 


---- 线性 Linear (灌溉 像 元 Irrigated pixels) 
0.60『 ee 线性 Linear (未 灌溉 像 元 Non-irrigated pixels) 
灌溉 线 Irrigated line Oo By 


NDYVI 
be 
全 
be 


3 ”基于 NDVI 与 NDWI 线性 关系 的 灌溉 线 确定 方法 
Fig.3 Method of irrigation line determination based on linear 
relationship between NDVI and NDWI 

NDWI: 归 一 化 差 值 水 分 指数 Normalized difference water index. 


1.5.4 ”精度 验证 
灌溉 像 元 提取 精度 ( a ) 采 用 下 式 计算 : 


N 
cw =(1- 一 人 )x1009%4 7 
( 和 0 (7) 


式 中 : NM 为 真实 灌溉 像 元 数 ; ,为 试验 中 提取 灌溉 像 
元 时 , 错误 分 到 未 灌溉 像 元 这 一 类 的 数量 错 分 数 。 

NDVI 校正 精度 验证 , 采取 A 地 块 试验 , B 地 块 
进行 验证 的 原则 ,通过 校正 精度 a 进行 判定 ， 且 对 上 比 
NDVI 校正 后 散 点 图 离散 程度 。 


2 结果 与 分 析 
2.1 裸 土 期 与 生长 期 灌溉 与 否 像 元 的 差异 比较 分 析 
4a 和 图 4b 分 别 显示 了 裸 十 期 (2002 年 5 月 30 
日 ) 和 棉田 生长 期 长 势 中 等 像 元 棉花 灌溉 与 未 灌溉 各 
波段 反射 率 差异 图 , 图 4c 为 灌溉 与 未 灌溉 像 元 4 种 
植被 指数 的 变化 差异 图 (2002 年 6 月 23 日 )。 从 图 4a 
可 以 看 出 ， 裸 土 时 期 各 波段 反射 率 差异 较 大 ,说明 
灌溉 会 造成 土壤 背景 发 生 显著 变化 ; 从 图 4b 可 以 看 
出 , 生长 期 中 等 长 势 像 元 的 灌溉 与 未 灌溉 像 元 各 波 
段 反射 率 变化 仍 较 显著 , 说 明 灌 溉 对 第 174 d 冠 层 棉 
花 反射 率 依然 存在 影响 ， 造 成 各 波段 反射 率 明 显 变 
化 。 从 4 种 VI 散 点 图 可 以 看 出 (图 4c), 4 种 VI 灌溉 
像 元 均 高 于 未 灌溉 像 元 , NDWI 变化 极 小 , 而 NDVI 
变化 较 大 , 约 为 12%， 其 他 3 种 土壤 调节 植被 指数 
变化 较 NDVI 更 大 ,没有 从 根本 上 剔除 土壤 背景 的 
影响 ， 不 能 真实 反映 植被 的 真实 长 势 。 


2.2 ”棉花 灌溉 与 未 灌溉 像 元 NDWI-VI 线性 回归 关 

系 分 析 

根据 灌溉 与 未 灌溉 像 元 的 4 种 VI 散 点 图 (图 4c) 
发 现 NDWI 变化 差异 极 小 ， 说 明 灌 溉 对 其 造成 影响 
很 小 。 因此 本 文 拟 以 NDWI 为 自 变 量 对 其 与 3 种 VI 
的 定量 关系 进行 研究 分 析 ( 图 5)。 从 图 5a、5b 可 以 
发 现 ， 灌溉 与 否 NDVI 与 MSAVTI 值 发 生变 化 , 即 灌 
溉 棉花 像 元 NDVI 与 MSAVTI 值 明显 高 于 未 灌溉 像 元 ， 
差异 在 10% 以 上 。 且 从 图 5 可 以 发 现 灌溉 与 否 棉花 
地 块 的 NDWI 与 3 种 VI 均 存在 极 显 著 的 线性 关系 。 
以 此 进行 线性 回归 分 析 ， 线 性 回归 模型 如 表 1。A、 
B 两 个 地 块 3 种 VI 与 NDWI 均 是 极 显著 相关 , A 的 
决定 系数 更 高 ,在 0.93 以 上 , B 的 决定 系数 为 0.90。 
由 此 可 以 看 出 ,灌溉 会 引起 植被 指数 变化 ， 且 变化 
差异 较 大 ,最 高 达 12.26%。 因 此 , 在 农作物 遥感 长 
势 精准 监测 中 , 保证 NDVI 值 实时 监测 的 准确 性 极 
2.3 ”基于 STDwi 方 法 和 IRL 方法 的 棉花 像 元 灌溉 

提取 及 精度 评价 

本 文 对 研究 区 种 植 棉花 的 A 地 块 进行 灌溉 像 
元 提取 ， 分 别 依据 IR_L( 灌 溉 线 法 ) 和 STDwi( 标 准 
差 方 法 ) 两 种 方法 ， 提 取 精 度 见 表 2。 IR_L 方法 精度 
可 达 94.72%， 相 比 STDwr 方 法 ， 更 适用 于 灌溉 信息 
提取 。 

本 文选 择 利用 A 地 块 试 验 , B 地 块 验证 的 精度 
验证 方法 。 再 次 选取 IR_L 方法 对 B 地 块 进行 验证 ， 
根据 结果 发 现 ， 本 文 提出 的 IR_L 灌溉 提取 方法 对 B 
地 块 依然 适用 ,精度 高 达 88.41%。 

2.4 作物 生长 初期 棉花 灌溉 信息 遥感 校正 

6 为 A 地 块 经 过 灌溉 校正 方法 调整 后 的 
NDWI 与 NDVI 的 散 点 图 , 与 未 消除 灌溉 影响 的 
NDWI-NDVI 散 点 图 (图 5a) 相 上 比较, 利用 IR_L 方法 
校正 后 的 NDVI 明 显 差异 减 小 。 且 经 过 校正 后 , NDVI 
差异 由 灌溉 前 的 12% 减 小 至 2%， 整 个 地 块 散 点 趋势 
趋 于 一 致 。 

3 讨论 

农作物 长 势 监测 主要 是 为 田间 管理 提供 及 时 准 
确 的 作物 生长 信息 ， 并 为 早期 估 产 提供 依据 中 。 作 
物 长 势 遥 感 监测 充分 体现 了 遥感 技术 全 面 、 及 时 、 
经 济 的 特点 。 近 十 几 年 , 许多 研究 利用 高 时 空 分 辨 
率 的 穆 感 数据 可 以 迅速 获取 大 范围 的 作物 长 势 信 
息 。 但 现实 情况 下 ， 由 于 受 大 气 和 土壤 背景 的 影响 ， 
传感器 只 收 到 来 自 目标 的 部 分 信号 ， 同 时 收 到 部 分 
噪音 , NDVI 也 由 此 受到 影响 ， 且 覆盖 度 越 低 ， 土壤 
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4 ”棉田 灌溉 与 未 灌流 像 元 的 反射 率 差异 图 
Fig.4 Reflectance differences between irrigated and 
non-irrigated pixels of cotton field 

a 裸 土 期 灌溉 与 未 灌溉 像 元 各 波段 反射 率 差 异 ; b: 生长 初期 ( 积 日 
174 d) 灌 溉 与 未 灌溉 像 元 各 波段 反射 率 差异 图 ; c: 生长 初期 ( 积 日 174 d) 
灌溉 与 未 灌溉 像 元 的 4 种 植被 指数 的 散 点 图 。NDWI: 归 一 化 差 值 水 分 
外 数 ; MSAVI: 调整 型 土壤 调节 植被 指数 ; SAVI: 土壤 调节 植被 指数 。a: 
band reflectance of irrigated pixel and non-irrigated pixel during bare soil 
period; b: band reflectance of irrigated pixel and non-irrigated pixel at early 
growing period (DOY 174); c: four vegetation indexes of irrigated pixel 
and non-irrigated pixel at early growing period (DOY 174). NDWI: nor- 
malized difference water index; MSAVI: modified soil adjusted vegeta- 
tion index; SAVI: soil adjusted vegetation index. 


Fig. 3 Two dimensional scatter diagrams of different vegeta- 
tion indexes and NDWI in the early growth period of cotton 
(DOY 174 d) 

a; A 地 块 NDVI-NDWTI; b: A 地 块 MSAVI-NDWI; c: B 地 块 
NDVI-NDWI。NDWI: 归 一 化 差 值 水 分 指数 ; MSAVI: 调整 型 土壤 调 
节 植 被 指数 。a: NDVI-NDWI scatter diagram of plot A; b: plot 


MSAVI-NDWI scatter diagram of plot A; c: NDVI-NDWI scatter dia- 
gram of plot B. NDWI: normalized difference water index; MSAVI: 
modified soil adjusted vegetation index. 


背景 对 NDVI 的 影响 越 大 。 本 文选 择 积 日 174 d 作 
为 研究 时 相 , NDVI 值 小 于 0.5, 仍 处 于 作物 生长 初 
期 ,由 NDVI 与 LAI 的 相关 性 可 知 覆 盖 度 较 低 时 , 土 
壤 背 景 对 NDVI 影响 较 大 ， 因 此 本 文 在 作物 生长 初 
期 对 由 于 灌溉 造成 的 土壤 背景 影响 进行 分 析 ， 对 
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表 1 灌溉 与 未 灌溉 像 元 的 归 一 化 差 值 水 分 指数 (NDWI, y) 与 植被 指数 的 线性 回归 模型 
Table ] Linear regression models of normalized difference water index (NDWI, y) and vegetation indexes for irrigated and 
non-irrigated pixels of cotton field plots 
Independent < 本 ER Vegetation index Pixels 
Plot variable 灌溉 Irrigation 未 灌溉 Non-irrigation change rate (%) All pixels ber 
y=1.236 2NDVI+0.291 2 y=1.226 7NDVI+0.259 4 y=1.372 1INDVI+0.258 2 $519 
A ND R=0.938 7 R’=0.973 7 1226 R2=0.945 9 
SAVI y=1.847 1SAVI+0.434 3 y=1.832 4SAVI+0.387 3 12.14 y=2.048 3SAVI+0.385 5 
R=0.938 8 R’=0.973 8 “ 及 2=0.9406 3 
J 一 1.178 SMSAVI+0.463 7 y=1.263 3MSAVI+0.421 2 y=1.35MSAVI+0.425 3 
MSAVI R=0.936 2 R=0.971 2 10:09 R=0.948 
y=1.353 8NDVI+0.254 6 y=1.403 2NDVI+0.216 3 y=1.370 8NDVI+0.238 4 
We R’=0.829 3 R’=0.889 2 EE R=0.753 90 


SAVI: 土壤 调节 植被 指数 ; MSAVI: 调整 型 土壤 调节 植被 指数 。SAVI: soil adjusted vegetation index; MSAVI: modified soil adjusted vege- 
tation index. 


表 2 基于 灌流 提取 方法 的 像 元 灌流 信息 提取 精度 


Table 2 Extraction accuracy of pixel irrigation information based on irrigation extraction method of cotton field plots 


方法 地 块 分 阶段 阅 值 灌溉 线 方程 灌溉 提取 精度 灌溉 像 元 数 总 像 元 数 
Method Plot Stage threshold Irrigation line equation Accuracy (%) Number of irrigation pixels All pixels 
STDwI A 0.029 48 85.09 322 519 
B 0.013 108 60.87 276 506 
IRL A y=1.231 SNDVI+0.275 3 17 94.72 322 519 
B y=1.378 SNDVI+0.228 7 32 88.41 276 506 


STDwi: 标准 差 法 ; IR_L: 灌溉 线 法 。STDwr standard deviation of NDWI (normalized difference water index) method; IR_L: irrigation line 


extraction method. 


0.60[ 口 灌溉 像 元 校正 后 Irrigated pixels after correction 


| 慷 节 参 数 过 的 确定 是 一 个 难题 ， 针 对 不 同 植被 覆盖 度 ， 
0.551 ^ 未 灌溉 像 元 Non-irrigated pixels 


L 值 的 确定 存在 难度 ， 植 被 覆盖 度 划分 并 不 明显 ; 
MSAVI 模型 中 没有 明确 反映 土壤 背景 的 参数 。 然 而 ， 
实际 上 的 植被 冠 层 光谱 是 一 个 组 合 ， 它 是 由 几 个 复 
杂 的 光学 物理 过 程 组 合 而 成 ， 仅 通过 线性 组 合 对 其 
进行 反 演 ， 存 在 一 定 的 局 限 性 。 因 此 很 难 在 实际 应 
用 中 通过 植被 指数 模型 定量 消除 土壤 背景 的 影响 。 

本 文 利用 变化 差异 最 小 的 NDWI 与 NDVI 的 线性 关 
系 ， 对 灌溉 后 的 NDVI 再 次 进行 校正 , 一 定 程度 上 
消除 低 覆 盖 度 下 土壤 背景 的 影响 ， 有 效 降 低 了 由 于 
灌溉 引起 的 土壤 背景 水 分 信息 的 影响 ,为 遥感 定量 
监测 提供 了 便利 ， 可 更 加 准确 地 对 作物 长 势 进行 预 
判 ， 推 进 精 准 灌溉 ， 节 约 水 资源 ， 促 进 农 业 可 持续 
发 展 。 

4 结论 


本 文 以 美国 加 州 棉花 为 研究 对 象 ， 对 研究 区 棉 
花生 长 初期 灌溉 与 未 灌溉 像 元 的 植被 指数 差异 与 各 
波段 反射 率 差 异 进行 分 析 ,， 利用 灌溉 与 未 灌溉 像 元 
的 NDVI 与 NDWI 之 间 极 显著 的 线性 关系 ,构建 两 
种 灌溉 提取 方法 ， 实 现 了 快速 准确 地 对 生长 初期 棉 
花 灌 溉 像 元 进行 提取 ， 且 精度 高 达 88%。 并 利用 灌 
溉 线 模 型 对 含有 生长 初期 棉花 灌溉 像 元 NDVI 值 进 
行 校正 ， 有 效 地 去 除了 土壤 水 分 对 植被 指数 NDVI 
的 影响 。 本 文 研 究 结 果 解 决 了 灌溉 造成 NDVI 值 
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0 0.05 0.10 0.15 0.20 0.25 
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6 棉花 A 地 块 NDVI 灌溉 提取 校正 后 散 点 图 
Fig.6 NDYVI scatter of Plot A of cotton field after NDVI irri- 
gation extraction and correction 

NDWI: 归 一 化 差 值 水 分 指数 。 NDWI: normalized difference 
water index. 


NDVI 进行 遥感 提取 与 校正 ， 这 对 作物 长 势 监测 有 
一 定 意 义 ， 同 时 有 利于 进一步 实现 精准 灌溉 。 

本 文 对 比 NDVI、NDWI、SAVI、MSAVI 4 种 
植被 指数 在 灌溉 前 后 各 波段 反射 率 的 差异 ,结果 显 
示 NDWI 在 灌溉 前 后 差异 最 小 ; 且 发 现 土壤 调节 植 
被 指数 系列 并 不 能 完全 消除 土壤 背景 的 影响 ， 尤 其 
针对 本 文 研 究 区 的 灌溉 情况 , 土壤 调节 植被 指数 消 
除 水 分 背景 影响 的 效果 并 不 好 , 与 NDVI 相 比 差异 
更 大 。 土 壤 调 节 植 被 指数 系列 多 数 是 以 土壤 线 模型 
计算 而 形成 的 线性 模型 ， 但 现实 中 土壤 线 难 以 获 
取 。 本 文选 取 的 土壤 调节 植被 指数 , SAVI 中 土壤 调 
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高 于 正常 水 平 的 问题 ,反映 了 植被 的 真实 信息 ， 
实现 了 对 生长 初期 作物 长 势 的 准确 遥感 监测 。 本 
文 研究 成 果 将 为 低 覆 盖 度 下 NDVI 的 研究 应 用 提 
供 参 考 依据 ， 也 将 会 促进 遥感 定量 监测 与 精准 农 
业 的 发 展 。 
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